Long-range ferromagnetic correlations between Anderson impurities in a 

semiconductor host 
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We study the two-impurity Anderson model for a semiconductor host using the quantum Monte 
Carlo technique. We find that the impurity spins exhibit ferromagnetic correlations with a range 
which can be much more enhanced than in a half-filled metallic band. In particular, the range is 
longest when the Fermi level is located above the top of the valence band and decreases as the im- 
purity bound state becomes occupied. Comparisons with the photoemission and optical absorption 
experiments suggest that this model captures the basic electronic structure of Cai^Mn^ As, the 
prototypical dilute magnetic semiconductor (DMS). These numerical results might also be useful 
for synthesizing DMS or dilute-oxide ferromagnets with higher Curie temperatures. 

PACS numbers: 75.50.Pp, 75.30.Hx, 75.40.Mg, 71.55.-i 
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The discovery of ferromagnetism in alloys of III-V 
semiconductors with Mn started an intense research ac- 
tivity in the field of dilute magnetic semiconductors 
(DMS) P, 0, El- Room-temperature ferromagnetism in 
DMS would be a significant development for spintronics 
device applications. In this respect, it is important to 
understand the nature of the correlations which develop 
between magnetic impurities in semiconductors and how 
they differ from that in a metallic host. With this pur- 
pose, we present exact numerical results on the two- 
impurity Anderson model for a semiconductor host. 

In order to study the multiple charge states of Au im- 
purities in Ge, the single-impurity Anderson model of 
a metallic host was extended to the case of a semicon- 
ductor host using the Hartree-Fock (HF) approximation 
After the discovery of DMS, the magneticproperties 
of this model were addressed within HF [H, [||. It was 
shown that long-range ferromagnetic (FM) correlations 
develop between Anderson impurities in a semiconduc- 
tor when the Fermi level is located between the top of 
the valence band and the impurity bound state (IBS), as 
illustrated in Fig. 1. The FM interaction between the 
impurities is mediated by the impurity-induced polariza- 
tion of the valence electron spins, which are antiferromag- 
netically coupled to the impurity moments. Within HF, 
the impurity-host hybridization also induces host split-off 
states at the same energy as the IBS. When the split-off 
state becomes occupied, the spin polarizations of the va- 
lence band and of the split-off state cancel. This causes 
the long-range FM correlations between the impurities 
to vanish. Within the context of DMS, the Anderson 
Hamiltonian for a semiconductor host was also consid- 
ered by Krstajic et al. [7|, and it was shown that an FM 
interaction is generated between the impurities due to 
kinematic exchange. In addition, this model was studied 
within HF for investigating the multiple charge and spin 
states of transition- metal atoms in hemoprotein Q. 

In this paper, we present quantum Monte Carlo 



(QMC) data on the two-impurity Anderson model for 
a semiconductor host and make comparisons with the 
HF results. We find that in a semiconductor the nature 
of the magnetic correlations between the impurities is 
different than in a metallic host. In particular, the impu- 
rities exhibit long-range FM correlations when the Fermi 
level is located above the top of the valence band, and 
the FM correlations weaken as the IBS becomes occu- 
pied in agreement with HF [|| . Comparisons with the 
photoemission and optical absorption experiments sug- 
gest that this model captures the basic electronic struc- 
ture of Gai^^Mn^As. These numerical results outline 
the parameter regime which yields the longest-range FM 
correlations, and this information might be useful for syn- 
thesizing DMS or dilute-oxide ferromagnets with higher 
Curie temperatures. 

The two-impurity Anderson model for a semiconductor 
host is defined by 
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where cj cct(T (ck QC r) creates (annihilates) a host electron 
with wavevector k and spin a in the valence (a — v) or 
conduction (a — c) band, d\ a (dj CT ) is the creation (anni- 
hilation) operator for a localized electron at impurity site 
i, and n^a = d\ a di a . The hybridization matrix element 
is Vuj = l^exp(ik • Rj), where Rj is the coordinate of 
the impurity site j. As usual, Ed is the d-level energy, 
U is the onsite Coulomb repulsion, and p the chemical 
potential. The valence and the conduction bands have 
the forms e£ = — D (k/k ) 2 and e£ = D (k/ko) 2 + Aq, re- 
spectively (Figure 1). Here, D is the bandwidth, ko is the 
maximum wavevector and Aq is the semiconductor gap. 
In this paper, we consider a two-dimensional semiconduc- 
tor host with a constant density of states po — k^/(ATiD). 



FIG. 1: (color online) Schematic drawing of the semiconduc- 
tor host bands ej? (solid curves) and the impurity bound states 
(thick arrows) obtained within HF. The dashed line denotes 
the chemical potential fx. 



We find similar results for the three-dimensional case . 
We determine the energy scale by setting D = 12.0. In 
addition, we use U = 4.0 and Ed = ix — U/2 (the sym- 
metric case) , so that the impurity sites develop large mo- 
ments. We note that this Hamiltonian is particle-hole 
symmetric with respect to half- filling at fx = Ag/2. For 
the DMS materials, it is estimated that Aq/D w 0.1 to 
0.2. We report results for Aq = 2.0, values of the hy- 
bridization parameter A = irpoV 2 ranging from 1.0 to 
4.0, and inverse temperature j3 = 1/T from 4 to 32. In 
order to study the evolution of the magnetic correlations 
as we go from a metallic to a semiconductor host, we will 
present results for fx from —D/2 to Ag/2. 

The numerical results presented here were obtained 
with the Hirsch-Fye quantum Monte Carlo technique 
[icj |. In the following, we will first present results on 
the impurity equal-time magnetic correlation function 



(MfM| 



where Mf 



n l( i] - n ldL is the impurity 



magnetization operator. Next, we will present results 
on the impurity single-particle spectral weight A(uj) = 
— (l/7r)ImGfj(w), which is obtained with the maximum- 
entropy analytic continuation technique [l3| from the 
QMC data on the impurity Green's function 



GUr) = -(T T d ta (r)dl(0)). 
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Here, T T is the Matsubara time-ordering operator and 
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Since the maximum-entropy pro- 



cedure requires QMC data with very good statistics, our 
results on A(co) will be limited to the high-temperature 
(3 = 8 case. For lower T, we will discuss QMC results on 
the impurity occupation number (rid) = (nidi) + ( n idi)- 
We will also show data on the zero-frequency inter- 
impurity magnetic susceptibility defined by 



Xi 2 (w = 0) 



dT(M 1 (r)M 2 (0)). 



(3) 
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FIG. 2: (color online) (MfMf) vs k R plotted at j3 = 16 and 
various fi for hybridization (a) A = 1.0, (b) 2.0, and (c) 4.0. 



The following results were obtained using Matsubara 
time steps At = 0.125 and 0.25. 

Figures 2(a)-(c) show the impurity magnetic correla- 
tion function (MfM|) vs k R, where R = |Ri - R 2 | is 
the impurity separation, at /3 = 16 for fi varying from 
-6.0 to 1.0. Fig. 2(a) is for hybridization A = 1.0. At 
fx = —6.0, we observe oscillations in the R dependence 
due to an RKKY-type effective interaction between the 
impurities. These results are similar to what has been 
obtained previously with QMC for a half-filled metallic 
band Till . HH ■ The wavelength of the oscillations in- 
creases when /j, moves to —1.0, because of the shortening 
of the Fermi wavevector. When fx = 0.0, the impurity 
spins exhibit long-range FM correlations at this temper- 
ature. We observe that upon further increasing \x to 0.5 
or 1.0, the FM correlations become weaker. This is be- 
cause the IBS becomes occupied as fx changes from 0.0 
to 0.5, as will be seen in Fig. 3(a). In Figs. 2(b) and 
(c), results on (MfM%) are shown for A = 2.0 and 4.0, 
respectively. In Fig. 2(b), we observe that (MfM|) has 
the slowest decay for [x = 0.5, while in Fig. 2(c) this 
occurs for [x = 1.0. We find that the impurity occupa- 
tion (rid) increases between fx = 0.5 and 1.0 for A = 2.0 
and (3 = 16. In addition, for A = 4.0 and f3 = 8, the 
maximum-entropy image of A(u>) shows that the IBS is 
located at w w 1.0. Hence, we observe that the range of 
the FM correlations for the semiconductor is determined 
by the occupation of the IBS in agreement with the HF 
predictions ]^,@]. In Figs. 2(a)-(c), it is also seen that 
the range increases with decreasing A. 

In Figs. 3 and 4, we discuss the A = 1.0 case in more 
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FIG. 3: (color online) (a) Impurity single-particle spectral 
weight A(lu) vs w for fco-R = 5 and 10 at /3 = 8. Here, the 
vertical dashed line denotes n, and the top of the valence band 
is located at ui = 0. In (b) and (c), (Mf Mf) vs p is plotted 
for fcoR = 5 and 10, respectively, at various (5. These results 
are for A = 1.0. 



detail. In Fig. 3(a), the impurity spectral weight A(uj) 
vs u> is plotted for j3 = 8, fj, = 0.1, and fco-R = 5 and 
10. Here, the w-axis has been shifted so that the top of 
the valence band is located at ui = 0. For fco-R = 10, we 
observe a peak at ujbs ~ 0.1 in the semiconductor gap, 
which we identify as the IBS. For fco-R — 5, the bound 
state is broader due to stronger correlations between the 
impurities. However, we also find that A(u>) exhibits sig- 
nificant T dependence at j3 = 8, and Fig. 3(a) does not 
yet represent the low-T limit. Next, in Figs. 3(b) and 
(c), (MfMf) evaluated at koR — 5 and 10 is plotted 
as a function of /i. Fig. 3(b) shows that, at low T for 
fco-R = 5, (MfMf) decreases when \i > 0.25. For this 
value of fco-R and (3 — 32, we find that the impurity oc- 
cupation (rid) develops a step discontinuity at fj, w 0.25, 
which is consistent with the decrease of (Mf M|) when 
ft > 0.25. For fc .R = 10 and (3 = 32, both (M{M%) 
and (rid) exhibit significant T dependence in the vicinity 
of the semiconductor gap edge. These results show that 
(Mf Mf ) depends strongly on the value of \i. 

Figure 4(a) shows the temperature dependence of 
(Mf Mf ) vs fco-R for fx = 0.1. We observe that, at = 32, 
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FIG. 4: (color online) (a) (MfMJ) vs fc R for /u = 0.1 at 
various (3. (b) Inter-impurity magnetic susceptibility Xi2(<^ = 
0) vs /3 for koR = 10 at various fi. These results are for 
A = 1.0. 



the range of the FM correlations is enhanced by about 
an order of magnitude with respect to that in a half- 
filled metallic band. In Fig. 4(b), the T dependence of 
the inter-impurity susceptibility Xi2( w — 0) is shown for 
fco-R = 10 and various values of fi. We observe that, as 
j3 increases, \i2 becomes strongly enhanced for /i = 0.1, 
while it remains weak for /j, = 0.3. At fi — —0.2, X12 is 
antifcrromagnetic for this value of fco-R. Figures 3 and 
4 show that long-range FM correlations develop between 
the impurities depending on the position of fi. 

In the QMC simulations, we find that (MfMf } has 
larger error bars when IBS is occupied. This might be 
due to the cancellation of the host spin polarizations 
originating from the split-off state and from the valence 
band. The intra-impurity measurements such as ((Mf ) 2 ) 
or Gfj(r) do not exhibit such behavior. 

Within HF, the range £ of the FM correlations be- 
tween the impurities is determined by the energy of 
the IBS and, hence, the energy of the split-off state: 
£ pa (16np uJBs)~ 1 ^ 2 for a constant density of states and 
semi-infinite host bands. In particular, the spatial extent 
of the spin polarization of the valence band around the 
impurity is given by 4poAZde~ r /^/(r/£), where Zd is the 
weight of the pole at cubs in the impurity single-particle 
Green's function. Within HF, ujbs decreases rapidly as 
A/Ag — ► 0, which agrees with the A dependence of the 
range seen in Figs. 2(a)-(c). In addition, for T = 0, 
A = 1.0, and < /Lt < uj B s, HF yields uj B s ~ 0.02 and 
fcoC ~ 9. We note that the effects of the inter- impurity 
correlations on the impurity single-particle Green's func- 
tion are neglected within this approximation. 
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In the QMC and HF calculations, the location of the 
Fermi level with respect to the IBS is important; the FM 
correlations weaken as the IBS becomes occupied. Pho- 
toemission and optical measurements on Gai^Mn^As 
provide evidence that the occupation of the Mn-induced 
impurity band is similarly important for the magnetic 
properties of this prototypical DMS ferromagnet. Pho- 
toemission experiments [141 ] observed an Mn-induced 
state above the valence band and right below the Fermi 
level in Gai-^Mn^As. Clearly, inverse photoemission ex- 
periments are required to detect the unoccupied portion 
of the Mn-induced impurity band. STM experiments also 
observed the impurity band in this compound [la |. Re- 
cent optical-absorption measurements [161 ] , which show a 
redshift of the mid-infrared peak with Mn doping, pro- 
vide evidence that the Fermi level is located in the Mn- 
induced impurity band in Gai_ K Mn x As. Furthermore, as 
the impurity band becomes less occupied with Mn dop- 
ing, the Curie temperature T c increases in annealed sam- 
ples, which is in agreement with the QMC and HF results. 
The comparisons of these experiments and the numeri- 
cal results suggest that the Anderson Hamiltonian for a 
semiconductor host provides a basic electronic model for 
the DMS ferromagnets. 

These numerical results also suggest that, in addition 
to Mn substitution, a possible way of increasing T c in 
Gai_ x Mnj;As is to decrease the impurity-band occupa- 
tion by changing the semiconductor host material or by 
using additional dopants. Alternative ways of enhancing 
the FM correlations in this model is provided by vary- 
ing the hybridization parameter A or the semiconductor 
gap Aq- The QMC simulations show that £ increases 
as A goes from 4.0 to 1.0. Within HF, £ can take very 
large values as A/Ag decreases. Recently, T c 's exceed- 
ing the room temperature have been reported in dilute 
oxides such as ZnO and Ti02 with transition metal impu- 
rities 17, HI]. At this point, it is important to determine 



whether the Anderson model of a magnetic impurity ap- 
plies to this case. Obviously, experiments probing the 
electronic state of the FM dilute oxides are necessary to 
answer this question. 

The two-impurity Anderson model for a semiconduc- 
tor host might be oversimplified for describing the ferro- 
magnetism of the DMS. Our calculations are for spin-1/2 
Anderson impurities, and we have neglected the multi- 
orbital structure of the spin-5/2 Mn impurities. Hence, 
the effects of the atomic Hund's rule couplings are not in- 
cluded. In addition, we neglect the long-range Coulomb 
repulsion between the impurity and the host electrons. 
Keeping these caveats in mind, it is interesting to note 
that the theoretical studies of high-T c DMS ferromag- 
netism could have preceded the experimental discovery, 
if the Hirsch-Fye algorithm had been applied to a semi- 
conductor host when it was developed twenty years ago. 

In summary, we have presented QMC results to show 
that long-range FM correlations develop between mag- 



netic impurities in semiconductors. In particular, the 
FM correlations have the longest range when the Fermi 
level is located above the top of the valence band, and 
they weaken as the IBS becomes occupied. Hence, the 
position of the Fermi level with respect to the IBS plays 
a crucial role in determining the range of the FM cor- 
relations in agreement with HF. Comparisons with the 
photoemission and optical absorption experiments sug- 
gest that the two-impurity Anderson model in a semi- 
conductor host captures the basic electronic structure of 
Gai^^Mn^As. The numerical results presented here out- 
line the parameter regime which yields the longest-range 
FM correlations, and this information might be useful for 
synthesizing higher- T c DMS materials. 
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